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An algebraic criterion that is sufficient to establish the existence of certain a priori estimates for 
the solution of first-order homogeneous linear characteristic problems is derived. Estimates of such 
kind ensure the stability of the solutions under small variations of the data. Characteristic problems 
that satisfy this criterion are, in a sense, manifestly well posed. 
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I. INTRODUCTION 

Any system of hyperbolic partial differential equations 
can be written in a peculiarly wicked form, namely: 
in characteristic form 0. In order to do this, a one- 
parameter family of characteristic surfaces is chosen as 
the level surfaces of a coordinate u, referred to as the 
null coordinate or retarded time. Instead of evolution 
from one time- level to the next, one obtains evolution 
from one characteristic slice to the next. However, char- 
acteristic surfaces are special: they are the only data 
surfaces for which the standard Cauchy problem can not 
be solved, because the differential operator of the system 
is internal to these surfaces, failing to evolve some data 
out of the surface. Naturally, the equations "degener- 
ate" , turning into a system where not all the equations 
can be integrated forward in retarded time, but there are 
"rules" that allow one to solve the problem in a hierar- 
chical manner, zig-zagging back and forth between the 
equations 2j. These hierarchical rules are the hallmark 
of characteristic evolution, making it significantly differ- 
ent from an initial value problem. 

It has been known for quite some time Q that the 
characteristic "Cauchy problem" - to obtain a unique so- 
lution from data given on a characteristic surface- can be 
solved so long as the data are split into two separate sets: 
some data are given on the initial null slice, and the rest 
of the data are given on another slice that must be trans- 
verse to the initial null slice. An issue that has attracted 
less attention is: if the data are perturbed slightly, under 
what circumstances is the variation of the solution un- 
der control? Equivalently, will almost-zero data evolve 
into a solution that is also close to zero? We set up a 
framework in which to address this question by defining 
certain types of estimates of the solution in terms of the 
free data, after Q- Subsequently we derive an algebraic 
criterion that is sufficient to determine whether the solu- 
tions satisfy such an a priori estimate, thus establishing 
their stability with respect to small variations of the free 
data. This kind of stability is of relevance to numeri- 
cal applications. A prominent instance of the use of the 
characteristic problem for numerical applications is that 



of the simulation of gravitational waves by numerically 
integratingthe characteristic formulation of the Einstein 
equations [j. 

As argued extensively in , characteristic problems for 
which the estimate can be established may be considered 
to be well posed in the sense that for each set of data 
the solution exists, is unique and depends continuously 
on the free data. In addition, characteristic problems 
that satisfy the algebraic criterion developed here can 
be thought of as manifestly well posed. Manifest well- 
posedness in the sense defined here is to characteristic 
problems as symmetric hyperbolicity 5] is to initial- value 
problems. 

Section [n] describes first-order linear characteristic 
problems after Duff 2J. The estimates of interest are 
defined in Section ITTT1 where the algebraic criterion is de- 
veloped as well. Concluding remarks are offered in Sec- 
tion Cvl 



II. HOMOGENEOUS LINEAR 
CHARACTERISTIC PROBLEMS IN 
CANONICAL FORM 

Consider a generic homogeneous hyperbolic system 
of linear partial differential equations for m functions 
v = {v a ) of n variables y a , which can be written in 
matrix form as follows 



dy a 



Dv = 0, 



(1) 



where summation over repeated indices is understood. A 
characteristic surface Af is a surface given by 4>{y a ) = 
such that 



det A 



dy a 



= 0. 



(2) 



Denote by m the multiplicity of this characteristic sur- 
face (so that the rank of A a d(j)/dy a is n — m). Suppose T 
given by ip(y a ) = is another surface intersecting Af at a 
submanifold of dimension n — 2, whose further properties 
are to be determined. We choose a suitable coordinate 
system (u, x, x l ) i — 1 ... n — 2 for R n adapted to these 
two surfaces; i.e., such that 
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(3a) 
(3b) 
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In these coordinates {Q reads 

B u d u v + B x d x v + B i div + Dv = 0, (4) 

with 

r)rh 



B x = A 

B i = A a 



dy a ' 

dy al 

dx 1 

dy a ' 



(5a) 
(5b) 
(5c) 



By J2J, there are to linearly independent left null vectors 
Zi v \ and also m linearly independent right null vectors 
zi v ) (with v = 1 . . . to) of the matrix B u , namely 



Z{ V )B U = 0, 



B u z 



0. 



(6a) 
(6b) 



We choose the right null vectors to be orthonormal in the 
sense that 



a a _ r 



(7) 



Multiplying (@J on the left with Z(„) we find that to of the 
equations in the system do not involve derivatives with 
respect to u: 

~z {v) B x d x v + z^B^iV + z {v) Dv = 0. (8) 

Our aim is now to find a convenient transformation of 
variables that takes advantage of this split of the equa- 
tions. We start by noticing that, using the to right null 
vectors as the first to legs of an orthonormal basis e' a of 
R m , we have a unitary transformation 



e'a = S a pep 



(9) 



from the trivial basis = {(1,0, . . . , 0), . . . , (0, . . . , 0, 1)} 
to the new orthonormal basis, with S a jSp-y — S a p and 
such that S ua = z^ a for v = 1 ... to. The components 
of v — v a e a in the new orthonormal basis are 



v' a = S a pVp. 



(10) 



In particular, the first to components are the scalar prod- 
ucts of v with the right null vectors z^ , which we denote 
by 

= z^ )a v a = v'^ for n = 1 ... to. (11) 

Multiplying (0J on the left with S, the system transforms 
into 



B' a d a v' + D'v' = 0. 



(12) 



From now on the index a refers to the characteristic co- 
ordinates, namely: a — u,x,i. The matrices have trans- 
formed according to B' a = SB a S t and D' = SDS\ 
and is the transpose of S. Because (B u S t ) av = 



B^pS v p = B^pZ^p — for v < to, the matrix B' 
has a Jordan form with all vanishing coefficients in the 
first to columns. This means that the it-derivatives of 
the to new variables w v are not involved, and conse- 
quently, the remaing variables v' v with v > to are the 
normal variables of the problem, which we denote by 
q = (g^), /i = to + 1 ... 7i. We have thus split the new 
fundamental variables into 



v' = Sv = (wi . . . w m ,q m+ i, . . . , q n ). 



(13) 



Inverting l|10(l we have v a = Sp a v'p, which can be used 
into © to obtain a set of equations in the transformed 
variables: 



B x S t d x v l + z {v) B i S t d i v' + S^DS^ 1 = 0. (14) 



We'd like for these equations to be solvable for the x- 
derivatives of all the variables w^, that is: the ones that 
do not evolve out of the initial characteristic surface. The 
first to terms in each equation for fixed v are 

z {u)a {B x S t ) atl d x w ll = z^ a B^pZ^)pd x Wf, (15) 

Thus the set of to equations l|14J) can be solved for the to 
variables d x w^ if and only if 



det {z{ v )aBl z^ )0 ) ^ 



(16) 



This is a restriction on the choice of ip(y a )- For this re- 
striction to hold it is sufficient, but not necessary, that 
the level surfaces of ip(y a ) be non-characteristic. In many 
applications, the level surfaces of tp are chosen to be time- 
like. For now on we assume that 1)16(1 holds. This allows 
us to interpret the to variables w v as the null variables 
of the problem. 

We have shown that under very weak conditions for 
the surface T, the most general characteristic problem 
takes the following form 



N u d u q + N x d x q + DPdi-uf ■ 
d x w + L x d x q + L%v' 



N°v' = 
hlV = 



(17a) 
(17b) 



Clearly the null variables w can be redefined by w = w + 
L x q so that none of the Eqs. ltT7b|l contains x— derivatives 
of the normal variables. Additionally, since N u is non- 
singular, we can choose normal variables q = N u q. In 
terms of these special choices of null and normal vari- 
ables, Eqs. I|17a|l and l(17b|l assume what is referred to as 
the canonical form: 

duq + ^d^+N^-did+N^ = (18a) 
d x w + L i d i v + L°v = (18b) 

where v = (w,q). We refer to (|18a|l as the evolution 
equations, and to (|18bl) as the hypersurface equations. 
For a unique solution to exist, one must prescribe the 
values of w on the surface x — and the values of q on the 
surface u = 0. The solution can then be constructed in a 
hierarchical manner. Since q is a known source for (|18b|) 
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at u = 0, then w can be found on the entire surface u = 0. 
Once u? is known at u = 0, it can be used as a given 
source for (|18a|) in order to find the values of the normal 
variables q on the next surface at u — du. These are 
then used into (|18bl) to obtain w on the surface u = du. 
And so forth. In fact, Duff proves a theorem of existence 
and uniqueness of the solution given the canonical form 
of the characteristic problem . 

As an example, consider the following equations for 
four unknowns v as functions of four variables x a = 
(t,x,y,z): 

d t vi = d x v 2 + d y v 3 + d z V4, (19a) 
d t v 2 = d x vi, (19b) 
d t v 3 = d y vi, (19c) 
d t v 4 = d z vi. (19d) 

These equations constitute a first-order version of the 
wave equation in three spatial dimensions (if we intepret 
the variables v a as the derivatives of a single function /). 
However, this first-order version of the wave equation has 
characteristic speeds of (rest) in addition to 1 (light). 
The level surfaces of <fi = t — x are null planes, so they are 
characteristic and intersect the surfaces of fixed value of 
x. We change coordinates (t, x, y, z) — > (u, x, y, z) with 



which implies that d t 
system (|19|1 turns into 



u = t — x 

— > d„, and d T 



d u v 1 +d u v 2 = d x v 2 H 

d u v 2 + d u v l = dxV 1 , 

d u v 3 = dyV 1 , 

d u v 4 = dzV 1 . 

We can read off the matrix B u : 



dyV" 



d r . - d .„ 



d z v 4 



110 
110 
10 
1 



(20) 
The 

(21a) 
(21b) 
(21c) 
(21d) 



(22) 



which is obviously singular of rank 3, so we have m — 1 
in this example. We expect only one null variable, and 
three normal variables for this problem in canonical form. 
B u has only one right null vector z = 2~ 1 / 2 (1, —1, 0, 0), 
and only one left null vector z = 2~ 1 / 2 (1, — 1, 0, 0), 
which coincides with z because B u is symmet- 
ric. An orthonormal basis for R 4 can be chosen as 
{2~ 1 /2( lj _ 1; o, o), 2-V2(l, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1)}. 
So the unitary transformation is 



S = 







1 1 

V2 V2 

10 

1 / 



(23) 



d z q4 



The null variable of the problem is w = (vi — V2)/V% and 
the normal variables are q u = {{v\ + V2)/V2,V3,V4), in 
terms of which the system (12111 takes the almost canonical 
form: 

1 
1 

d u qi - ^^92 
d x w T^dyQa - 



2d u q 2 - d x q 2 - 
d u q3 



- : d y q 2 



1 



0, 



0, 



1 

V2 
1 



d z w = 0, 



V2 



V2 



d z q4 



0. 



(24a) 
(24b) 
(24c) 
(24d) 



For a unique solution, we need to prescribe the value 
of w on the surface x = 0, and the values of 92,93 and 
94 on the surface u = 0. Notice that, in this example, 
the surface x — is not timelikc with respect to the 
hyperbolic operator A a , but is also characteristic, as can 
be seen by inspection of the matrix B x : 



B a 



0-100 
-1 





(25) 



However, we have zB x z = 2^0. Therefore the condi- 
tion l|16fl is satisfied in spite of the fact that the surfaces 
of fixed value of x are not timelike. 



III. WELL-POSEDNESS OF HOMOGENEOUS 
LINEAR CHARACTERISTIC PROBLEMS 

The canonical system (I18fl can be written in the com- 
pact form 



C a d a v + Dv = 



(26) 



where C and C x have block-diagonal forms of a special 
type: 



C u 



1 




N x 
1 



(27) 



where 1 is the identity of dimension n — m in the case of 
C u , and of dimension m in the case of C x . The matrix 
N x is square, of dimension n — m, and the various rectan- 
gular blocks are vanishing matrices whose dimensions 
are clear from the context. Dropping the hats ("") for 
ease of notation, the variable v represents the set of nor- 
mal variables q and null variables w of the characteristic 
problem in canonical form. Additionally, m functions 
wo = (wq(u, x 1 )) are given as data on the surface T and 
n — m functions go = {lo( x i x *)) are gi ven as data on the 
surface N. 

For the remainder of this section, we make the strong 
assumption that N x and C % are symmetric. Multiplica- 
tion of l|26l) by v on the left then leads to a "conservation 
law" of the form 



d a {vC a v)+vRv = 



(28) 



4 



where R = 2D — d a C a . We now integrate this conserva- 
tion law in an appropriate volume V of 1". Our volume 
is a "hyperprism" limited by the surface u — from "be- 
low", the surface x — on the "left", and the surface 
u + x = T, for an arbitrary constant T, on the "top" . 
We assume there are no boundaries in the remaining co- 
ordinate directions, in the sense that the solutions v will 
either be periodic functions of x l or will decay sufficiently 
fast at large values of x % in order for the integrals of their 
squares to exist. The integration yields 

f v{C u + C x )v dzZ T - I {vC u v)dAf- 

- [ (vC x v) dT+ [ vRv dV = 0. (29a) 
Jt Jv 

Clearly 

/ vC u vdN=( V {q v fdAf = |M| 2 , (30) 



if— m+1 



and 



j vC x vdT = J qN x qdT + J J2«) 2dT 



T Jl 

qN x q dT + \\w 
r 



2 



(31a) 



Thus Eq. (|2*9"|) is rearranged to read 
/ v(C u + C x )v d^T = ||9o|| 2 + lko|| 2 + [ qN x qdT 

vRvdV. (32a) 



If N x is non-positive definite but also such that 1 + N x 
is positive definite, we can define the norm of the solution 
v on the surface by 



HI 



f v(C u + C x )v dYiT: 



(33) 



and Eq. (|^2l implies 



\v\\t < \\qo\\ 2 + 



\w 1 1 2 - / vRvdV. (34) 
Jv 



In special case of constant coefficients with no undiffer- 
entiated terms, namely R = 0, Eq. I|34|) takes the form 



INt < IMI 2 + IKIf 



(35) 



which represents an a priori estimate of the solution 
in terms of the free data. It implies that the "size" of 
the solution is controlled by the "size" of the data. We 
may interpret it as a statement of well-posedness of the 
characteristic problem. Clearly the estimate holds in the 



presence of non-constant coefficients and undifferentiated 
terms as long as R is non-negative definite. 

An estimate can still be drawn in the presence of a 
negative definite bounded R, but it is weaker and holds 
only for small values of T, as we show next. 

Since R is negative definite, then 



-vRv <r^2(v u ) 2 



(36) 



where r =max(|i?y |) in the volume V, assuming that 
such a number r exists. On the other hand, since 
C + C x is positive definite and symmetric then all its 
eigenvalues are positive and we have 



v(C u + C x )v > cJ2( v "f 



(37) 



i/=i 



with c being the smallest eigenvalue of C u + C x . This 
implies 



Thus 



-vRv < -v(C u + C x )v 

c 



(38) 



vRvdV < - / \\v\\f dt (39) 
'v c Jo 

where |H| 2 is the norm of the solution on the surface S t 
given by u + x = t for fixed value of t < T. Thus the 
inequality 1)34(1 implies 



MIt < lkolr + IKH + 



HI?*- (40) 



Here ||<?o|| 2 and ||wo|| 2 are the norms of the normal and 
null variables with respect to the surfaces Af and T both 
bounded by the spatial surface at u + x = T . For any 
value of t < T we can write down the same inequality 

IMI 2 </ E(9")Wt+ / twvWldt' 

J N't Jr t c J 

(41) 

where Aft and T t are the subsets of and T bounded 
by S t , respectively. Since both integrals indicated are 
less than the norms \\qo\\ 2 and ||u>o|| 2 respectively, this 
implies 

rt 



Ml?<ii9dir + iKir + 



M\l dt' (42) 



Using this inequality recursively into the right-hand side 
of (|40|) we have 



\V\\ 2 T < 1 



rT (rT) 2 
~c~ + 2c 2 



+ . 



{rry 



x( ||<7o|| 2 + lko|| 2 + 



+ (r/ C y +1 / dh / dt 2 ... \\v\\ 2 t . +1 dt j+1 



(43) 
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for any given non- negative integer j. In the limit for 
j : —>■ oo the sequence in the right-hand side converges if 
(rT/c) < 1, in which case we have 



MIt 



(lMI 2 + IKH 2 ) (44) 



This is our final estimate for the solution in terms of the 
free data on the surfaces N and T. The estimate involves 
an exponential factor essentially due to the presence of 
undifferentiated terms. The exponential factor depends 
on the properties of the system of equations (the princi- 
pal matrices and the undifferentiated terms), but not on 
the choice of data. This is analogous to the a priori esti- 
mates for Cauchy problems with undifferentiated terms. 
As usual in such cases, the estimate is useless for large T, 
and, in particular, our proof only guarantees the estimate 
for T < c/r. Perhaps with greater care the estimate can 
be extended to longer values of T. 

Because the a priori estimates (|44H are independent 
of the choice of data, we can say that our characteris- 
tic problem is well-posed. The conditions under which 
we are able to derive a priori estimates thus become our 
criteria for well-posedness of linear homogeneous charac- 
teristic problems in canonical form: 

i) The principal matrices C a are symmetric. 

ii) The normal block of the principal x— matrix, denoted 
N x , is non-positive definite but such that 1 + N x is 
positive definite (namely, — 1 < N x < 0). 

There is, clearly, no obstacle in generalizing the con- 
struction slightly to the case where the characteristic 
problem is cast into "almost" canonical form, namely, 
the case when the principal matrices are 



C u 



N u 




N x 
1 



(45) 



which corresponds to a strictly canonical form up to a 
transformation of normal variables among themselves. 
In this case, the criterion is 

i) The principal matrices C a are symmetric and 

ii) The normal block of the principal u— matrix, denoted 
N u , is positive definite. The normal block of the 
principal x— matrix, denoted N x , is non-positive definite 
but such that N u + N x is positive definite (namely, 
—N u < N x <0). 

If i) and ii) hold for a linear homogeneous characteristic 
problem in "almost" canonical form, then the problem is 
well posed in the sense that there exist a priori estimates 
of the kind < e KT (||g || 2 + IKI| 2 ), where K is 

a constant independent of the data. This inequality is 
sufficient to establish the stability of the solutions under 
small variations of the data. Notice that —N u < N x < 
is equivalent to the requirement that the surfaces given 
by (t>{y a ) + ?P(y a ) = T with fixed value of T are spatial 
with respect to the hyperbolic operator A a , which in turn 
means that they can be interpreted as the level surfaces 
of a time function t(y a ) = <fi(y a ) + ip(y a )- 



As an illustration, we can see that the first-order form 
of the wave equation, Eqs. (|24|l . is well posed. For 
Eqs. (|21J) we have 



AT = 



and also 




N x = 




(46) 



C y 




/ 1 
J_ 

•v/2 I 10 1 
\0 1 
(47) 

Therefore all the conditions are satisfied, and the esti- 
mates follow. In this case, the estimates are of the form 
(|35|l and hold for any chosen T because all the princi- 
pal matrices have constant coefficients and there are no 
undifferentiated terms (R = 0). 



IV. CONCLUDING REMARKS 

Characteristic problems for hyperbolic equations are 
rarely discussed in the literature. In fact, prior to 
Balean's work 0, 0, Q> practically nothing was known 
about the characteristic problem of the simplest hy- 
perbolic equation, that is, the wave equation. Balean 
discussed how to derive estimates for the solutions of 
the wave equation in its standard second-order form. 
Balean's estimates differ markedly from ours. The es- 
timates for general linear characteristic problems of the 
first-order that we present here constitute a direct gener- 
alization of the estimates that we recently derived for the 
particular case of solutions of the characteristic problem 
of the wave equation as a first-order system of PDE's . 

The value of the generalization that we present here 
resides in the formulation of algebraic criteria sufficient 
for the existence of the a priori estimates. We demon- 
strate elsewhere that these criteria allow us to formu- 
late the characteristic problem of the linearized Einstein 
equations in a form that is guaranteed to be well posed. 

Several issues of interest remain wide open. First, 
given a general characteristic problem that is well posed 
in the sense that we introduce here, it is not at all clear 
as yet whether estimates of the derivatives of the solution 
in terms of the derivatives of the data would exist as well. 
We have succeeded in deriving estimates for the deriva- 
tives in the particular case of the characteristic problem 
of the wave equation However, the derivation de- 
pends strongly on the particular form of the hyperbolic 
operator of the wave equation, and its generalization to 
arbitrary characteristic problems is far from straightfor- 
ward, quite unfortunately. 

Secondly, a sufficient criterion to establish well posed- 
ness of a characteristic problem is useful, but a necessary 
criterion would, perhaps, be invaluable as a means to rule 
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out unstable problems with an eye towards numerical ap- 
plications. 

Thirdly, whether or not all well-posed hyperbolic prob- 
lems admit well-posed characteristic problems in our 
sense might well be the most intriguing open question 
at this time. 
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